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Abstract. This paper purposes to study quasi-normal modes due to massive scalar 
fields. We, in particular, investigate the dependence of QNM frequencies on the field 
mass. By this research, we find that there are quasi-normal modes with arbitrarily long 
life when the field mass has special values. It is also found that QNM can disappear 
when the field mass exceed these values. 

> 

o 

O ■ PACS numbers: 04.20-q, 04.70-s 

o . 

^ ■ 1. Introduction 

qh| Quasi-normal mode (QNM) is one of the important and exciting themes in the black 
jli ! hole physics. QNMs represent the behaviours of fields on a black hole spacetime under 
a certain condition, imposed both at the event horizon of the black hole and at the 
asymptotic infinity of the spacetime. This means that QNMs can be strongly affected 
^ i by the curvature of the spacetime. In other words, we can expect that the property of 
5t ! QNMs strongly reflects the property of the spacetime, and so we can grasp the property 
of the spacetime by studying on QNMs. 

The existence of QNMs was mentioned first by Vishveshwara[lJ in the context 
of scattering problem by black holes. Later, people began to perceive that QNMs 
were important sources of gravitational waves, and many researchers developed various 
methods to calculate QNM frequencies j2], 0, 0, [O], IZ], |H]; Recently, some 
researchers have pointed out that the asymptotic behaviour of the distribution of QNMs 
relates to the quantum gravitational theory [TOj, [H]. Thus, we recognise that researches 
on properties of QNMs will further increase its importance. 

Our motivation of this article is to study properties of a black hole spacetime using 
QNMs as a mathematical tool. Almost all studies on QNMs discuss only massless 
fields such as electromagnetic waves and gravitational waves, since these studies mainly 
focuses on phenomenological significances. Therefore, there are few studies on QNMs 
by massive fields (hereafter, we call them massive QNMs). One of the studies related 
to massive QNMs is a work by Simone and Will ^2] , which investigates massive QNMs 
on black hole spacetime using WKB method. They study the dependency of QNM 
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frequencies on the mass of the field, but their discussion is restricted to narrow range of 
the field mass due to the restriction required by the WKB method. Hence, their analysis 
cannot fully reveal the dependency of QNM frequencies on the field mass. Another work 
on massive QNM is a study by Konoplya P^. He study massive QNMs on a charged 
black hole back ground and consider QNMs due to massive and charged fields. His paper 
mainly focuses on correlation between QNMs and the charge of the field and the black 
hole. Therefore, the dependency of QNMs on the field mass is not argued , though his 
article investigate almost same situation that we will study in this article. 

Hence, we will study behaviours of QNMs to the field mass in detail. For this 
purpose, we first consider a model in which we can derive QNM frequencies analytically. 
Through this investigation we find that there is a singular phenomenon that QNMs may 
disappear when the field mass becomes sufficiently large. Second, we study QNMs on 
the Reisner- Nordstrom black hole spacetime by a numerical method. As a result, we 
confirmed that the singular phenomenon observed in the model can also be observed in 
the real black hole case. 

This article consists of the following contents. In section |21 we study a correlation 
between QNM frequencies and the field mass in a model which can be solved analytically. 
In section El we examine the real black hole case with a numerical method. In 
section HI we discuss the result obtained in section El In section we summarise our 
investigation. Additionally, [Appendix A shows the review of the calculation algorithm 
used in section |3] 

2. Analytical model 

In this section, we study the dependence of massive QNM frequencies on the field mass 
by using an analytically solvable model. 

Generally, in order for QNMs to exist, it is necessary that the background has a 
potential that can trap waves. Then, the trapped waves will exude out gradually. This 
phenomenon is exactly QNMs. If the asymptotic values of the potential are the same at 
both side of the potential peak, the waves will be the massless QNMs. If the asymptotic 
values are different, they will be massive QNMs. 

Let us consider the following model: 



+ - V{x.) 
dxi 



Zix^ = Q ila) 



« = ^(l + tanh|^). (Ic) 

Here, the parameters in the above equations cannot have any physical meanings, but we 
can consider some correspondences to the black hole parameters. M and m correspond 
to the black hole mass and field mass, respectively. L determines the width of the 
peak and this is dependent on M in the black hole case. The coordinate variable 
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Figure 2. Profiles of tlie potential in a Schwarzschild black hole case (M = 1). 
It is normalised by m^. 



corresponds to the so-called "tortoise" coordinates. The potential of this model (fT^ is 
shown in figure ^ For comparison with a real black hole one, we also show the potential 
for a Schwarzschild black hole. Because both potentials are very similar, we can also 
expect that the behaviours of QNMs will resemble in both cases. 

This model can be solved analytically as follows. First, by changing the coordinate 
X* to M, (|Tajl can be rewritten into the following form. 



ul 



, d 

au u{l 



u 



Z{u) = 0. (2) 



4M2 l-u_ 

Note that the infinity of x^, = ±oo corresponds to m = 1 (+oo) and {—oo). These two 
points u = and 1 are the regular singularities of the differential equation (j2)). Then, 
the behaviour of Z{u) around these two singularities is given by 

(3) 
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where = — L^z/^, 6^ = — L^(z/^ — m^). Therefore, by setting 

Ziu) = u^il-ufgiu) (4) 
and substituting (jlj) into Q, we obtain a differential equation for g{u) as follows, 



(P d 
u{l - u)— + {(1 + 2a) - 2(1 + a + h)u} — 



2 



+ <! -a - 6 - 2a6 - L^rri^ - + 2^^^ 

4M^ 



^7(«) = 0. (5) 



This is exactly the hyper-geometric equation and the solution of this equation can be 
given by the hyper-geometric function F(a, /5, 7; u). In consequence, we find the solution 
of as follows, 

Z{u) = u''{l-ufF{a,f3,r.u), (6a) 
1 . , / f 



7 = 1 + 2a. (6c) 
Now, let us give our attention to the following relation 

= W^^^^r^y^'^ - «)'^("' /5, 1 + « + /? - 7; 1 - n) 

r(7 - a)r(7 - /?) 

+ ^^Vr ZL - ")^^(^ _«,7-/3,l-a-/5 + 7;l - u).{y) 

According to (|7j), we can analytically determine QNM frequencies from the behaviour of 
the solution around u = and 1, because QNM is defined by the boundary conditions 
at = ±00, i.e. u = 0, 1. The boundary condition is that: QNMs are 

(i) ingoing waves at x* ~ —00 (u = 0), 

(ii) outgoing waves at x* ~ +00 {u = 1), and 

(iii) damping modes %. 

From the first condition at m = 0, we have to choose the sign of a as, 

a = —iLv. (8) 

The second condition at m = 1 means 
r(7)r(a + /3-7) _ „ 

r(a)r(/5) 

because (1 — m)~^ represents the incoming waves at m = 1 when we set h = +iL\/ — 
§. Therefore, QNM frequencies can be determined by, 

a, /? =-n{n = 0,1,---), (10) 

X We suppose the time dependency is exp(— wi) 
§ The branch of is — 7r/2 < argy^ < 7r/2. 



Massive QNM 



5 



according to ©• By solving this equation for u, we obtain 

, (2n + 1) (4n^ + 4n + /i^ - Ajj^^e^) 
4L (4n2 + 4n + /x^) 



y-±ib -r '±/(, -r ^ -r c ; / . /- ^-x 



where we set L = /xM and e = mM. When /i < 1, z/ is pure imaginary. However, i/ 
becomes a complex number if > 1. In this case, the first term in the right hand 
side of flTT|l represents the damping part and the second term does the oscillation 
part. Consequently, the damping coefficient of QNMs approaches to zero as a quadratic 
function of e. However, u never becomes a real number because the branch of square 
root changes discontinuously when QNM frequencies get across the real axis. According 
to this result, we can conclude that there exist QNMs with arbitrarily long decaying 
time when m changes appropriately. In other words, we find an interesting fact that 
there is a kind of resonance mode in a limiting situation and QNMs can disappear when 
m exceeds a certain value. 

This disappearance is due to the relation of the height of the peak of the potential 
with the field mass. Therefore, we can expect that the same phenomenon will occur 
in the black hole case. Thus, we will investigate the massive QNMs on the black hole 
spacetime in the next section. 

3. Massive QNM on black hole spacetime 

In this section, we consider QNMs of a massive scalar field on a spherical black hole 
spacetime. In the black hole case, we must tackle the problem in a numerical way 
because it is difficult to study it fully analytically. 

The metric of a spherical symmetric spacetime is given by 

ds^ = -f{r)de + f-\r)dP + r^dn\ (12) 

In this article, we consider the Reisner-Nordstrom black hole as the back ground 
spacetime, so the metric function /(r) is given by: 

/x =1 + ^ = l-- + \, 13 

where M is the black hole mass and Q is the black hole charge. We also introduce the 
normalised variables x and q as below, 

X = ^. (14a) 
, = |. (14*) 

respectively. 

In this spacetime, the radial equation for a scalar field in term of the normalised 
variables is given by 



(P Id / z/2 /(/ + 1) m2 ^ 



dx 



dx"^ f dx \ P x'^f f xf 



Z{x) = 0, (15) 
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Figure 4. Potential in case of e = 0.3 



where u and e are normalised frequency and mass of the scalar field and they are related 
to original frequency cu and mass m || by 

u = Mlo, (16a) 
e = mM, {16b) 

respectively. 

In order to find the potential for the scalar field, let us introduce a so-called tortoise 
coordinate by 

dX^: 1 

dx 7' 

and we obtain the following equation, 



;i7) 



0, 



(18) 



The original field equation is (□ — m )(f) = 0, and u is introduced by variable separation, 



-iiot 
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where the potential V{x^) is given by 

V{x.) = f + + I ) X fix). (19) 

Note that x in the right hand side of (jl9|) is a function of x^ by (jl7|) . The profiles of 
this potential in some cases are shown in figure El and figure |3] Figure |21 in section El 
also shows the potential for comparison with the potential of the analytic model. 

In order to calculate the QNM frequencies, we solve ()15|) numerically. We consider 
the case of |g| < 1 in this article, because the mathematical structure of (fT3|) with q = 1 
is apparently different from the case of |g| < 1. The reason of excluding the case of 
\q\ > 1 is that there is no horizon in this case and then the central singularity is naked. 
In this case, QNMs cannot be defined. Although the Schwarzschild case, q = 0, seems 
to be different case, we verify that there is no difference between the two cases in our 
results: (i) g = and (ii) q 0. 

The numerical method which we use in calculating QNM frequencies is the 
continued fraction method developed by Leaver 0. In jAppendix A[ we will briefly review 
this method for our case. In this article, we examine only / = and calculate QNM 
frequencies from e = to around 0.4 in several cases of g = (Schwarzschild case), 0.25, 
0.5, 0.75, 0.95 and 0.99. The results are shown in figureEl figurelHl figuredand figurelHl 
In these figures, i>i and i^r mean the imaginary and real part of the QNM frequencies, 
respectively. 

Figures El and El show the dependence of the QNM frequencies on the mass of the 
scalar field. It is obvious from the figures that the imaginary part of the QNMs on 
the black hole spacetime approaches to as well as the analytical model in the last 
section. This tendency is more clear in figure [71 Figure [71 shows the trajectories of the 
QNM frequencies on the u-plane when the field mass e changes while keeping q constant. 
From this figure, we can find that the trajectories of the QNM frequencies are almost 
linear, and we can guess the frequencies when QNMs disappear. We call the modes 
with these real frequencies as quasi-resonance modes ( QRM). According to the linearity 
of the QNM trajectories, we can calculate the dependence of QRM frequencies on the 
field mass e with the following fitting function, 

^1 = k{^R- ^qrm)- (20) 

The behaviour of this z/qrm to e is plotted in figure El 

We will give more detailed discussion on our results in the next section. 

4. Discussion 

In this article, we have studied massive QNMs for the two cases, the analytically solvable 
model and the black hole. We can find the same phenomenon, that the life of QNMs 
becomes longer as the field mass increases, in both cases. From this fact, we can conclude 
that there exists a QNM, oscillating without damping for an arbitrary period. This 
extreme mode may be regarded as a kind of resonance mode (after this fact, we named 




Figure 5. Dependence of Im(j/) to the field mass e for tlie cases oi q — 0.0, 0.25, 0.50 
and 0.75. 
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Figure 6. Dependence of \a\[y) to the field mass e for the cases oi q — 0.95 and 0.99. 



this mode quasi-resonance mode, QRM, in the last section). Now, we propose a physical 
picture of the disappearance of QNM and the existence of QRMs. 

First, let us discuss the case of the analytic model. We consider only n = case 
for simplicity. By setting n = in we obtain 

T — (,VT^). (21) 

What is noticed from is that QNM frequencies are almost determined by the width 
L of the potential peak when e ~ 0, and they turn out to depend on as e increases. 
Since the energy of QNMs is proportional to the real part of z/, the larger the energy 
of QNMs is, the gentler the peak of the potential becomes. Then the peak will vanish 
at e = 1/2. According to the WKB analysis of QNM by Schutz and Will g], QNMs 
can be regarded as the waves trapped by the peak of the potential. Therefore, in their 
picture, vanishing of the peak is equivalent to the potential being unable to trap any 
waves. This means that QNMs disappear. 
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Figure 7. Trajectory of QNM frequencies in z/-plane. 
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Figure 8. Dependence of QRM frequencies on the charge of the black hole. 
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Figure 9. Fitting by a parabola = cq + C2e^. The solid line shows the QNM 
trajectory of q = and the dashing line shows the fitting curve. 




Massive QNM 



10 



Next, let us consider the black hole case in this picture. In the black hole case, the 
peak of the potential vanishes at, for example, e = 0.25 for q = and e ~ 0.256 
for q = 0.5. On the other hand, figure El and figure El indicate disappearance of 
QNMs approximately at e = 0.4 ~ 0.5. Thus, the above picture seems to be realized 
approximately. Because the QNM frequencies are complex numbers, it is unreasonable 
to think that the disappearance of QNMs perfectly coincides with vanishing of the peak 
of the real potential in general. Hence, it is natural to consider that the coincidence 
in the analytic model is accidental. However, the behaviours of QNMs to m in both 
cases resemble each other very well. Figure El shows that the trajectory of QNMs can 
be fitted by a parabola very well. That is, our analytic model imitates the black hole 
case precisely. Therefore, since the analytic model predicts the disappearance of QNMs, 
it is undeniable that there exist QRMs in the case of black hole as well as the case of 
analytic model, and that QNMs will disappear as the field mass exceed a certain value. 

We also examine the effect of the black hole charge. We can see the behaviours of 
QNMs and QRMs when we change the black hole charge, q, in figure U\ and figure El 
These figures indicate that QNM and QRM frequencies tend to get larger when q 
increases. This can be explained by the potential profiles. According to figure El and 01 
when q increases, the peak of the potential becomes higher and the value of e at vanishing 
point of the peak becomes larger. Consequently, QNMs can exist at larger e and QRM 
frequencies get higher. These facts, too, support our picture of the disappearance of 
QNMs. 

5. Summary 

Let us summarise our results. First, we studied massive QNMs in the analytical model, 
and derived QNM frequencies analytically. By this analysis, we found that QNMs will 
disappear when the field mass increases and exceeds a certain value. Secondly, we 
examined the dependence of QNM frequencies on the field mass by numerical method 
in the case of black hole. There, we again found the disappearance of QNMs as we did 
in the analytic model case. We also studied the dependence of QNM frequencies on the 
charge of the black hole. Based on these examinations, we proposed a picture about the 
disappearance of QNMs. 

Appendix A. Continued fraction method 

In the appendix, we briefly review the continued fraction method, which is first 
developed by Leaver [3j. 

First, let x± be the solution of /(x) = in (fT5|l . Then, the points x = 0, X — , are 
regular singularities and the point x = oo is an irregular singularity of the differential 
equation ()15|) . Note also that x = x_ corresponds to the inner horizon of the black 
hole and x = x+ corresponds to the outer horizon. QNMs on a Reisner-Nordstrom 
black hole spacetime exist in regions connected to the spatial infinity of the spacetime. 
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Figure Al. Conformal diagram for a Reisner-Nordstrom black liole witli \q\ < 1. 



These regions are shown as A, A' in the conformal diagram ffigure IXT| . Therefore, the 
boundary conditions for QNMs on this spacetime are imposed ai x = x+ and at x = oo. 



Let us expand the solution of (|THjl as follows, 



Z{u) = e'^^'uPix - J2 

n=l 



where 



u 



2 2 2 

=X - e , 
x+^ 

P = -^x ^ _^ , 

6= -l + ^(z/2 + x')(x+ + a;_). 
Substituting (jA.lj) into (fTHj) . we obtain the following recursion relation, 
A\a\ + -Biflo = 0, 

A„a„ + i?„a„_i + C„a„_2 = 0, for n = 2, 3, 

where 

An = n^i^x^ — + 2mx+^x, 



'5 "5 ) 



Br, = — 2n (x„ — 



+ n X < — X-j 



[x- + x+) — iy(x- — 3s I 



+ 2(x- — x+ — 2ix+^x) 
+ <!-[/(/ + + x+) 



(A.i; 



(A.2) 
(A.3) 

(A.4) 
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+ [z/^(x_ - x+) + 2ix - (x_ + 3x+)x^]x+" 

,2 



- -{x- - x+ - 2ix+'^x) 



X 



[X- + 



uix- — SXj 



,(A.5) 



+ n X —ilu-\- 



r 



V 



Xj^) — 2(a;„ — x+ — ix+^x) 



j (a;_^ — a;+^) (a;_ + a;+) 



+ 



[x^ - x^) - Hxj^^x + 1^^- + + 3x+^)x^ f-'lA-e) 



Using these coefficients y4„, i?„ and C^, QNM frequencies are given by the vanishing 
point of the following continued fraction equation, 

AxC2 A2C3 An-iCn 



= Bi 



(A.7) 



B2— -B3— Bn — 

For practical calculation, we must choose n for the needed accuracy and solve the 
algebraic equation of finite but very large degrees. 
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